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Проведено исследование импедансных характеристик последовательного резонансного контура цепи нагрузки системы индук- 
ционного нагрева. Предложен алгоритм определения границ изменения резонансной частоты и добротности колебательного 
контура при изменении температуры нагреваемого объекта, величины зазора между индуктором и объектом, числа витков и то- 
ка индуктора, номинальной частоты. 


Введение 

В настоящее время высокочастотный индук- 
ционный нагрев активно применяется в промы- 
шленности для широкого спектра технологических 
операций. Источником магнитного поля, вызы- 
вающего нагрев детали, является индуктор. Обра- 
зованная в результате электромагнитного взаимо- 
действия с деталью система «индуктор - нагревае- 
мый объект» является для преобразователя частоты 
комплексной активно-индуктивной нагрузкой, 
индуктивный характер которой компенсируется 
конденсатором, с образованием последовательного 
резонансного контура. Важной особенностью си- 
стем индукционного нагрева является зависимость 
электрофизических свойств детали от ее темпера- 
туры, поэтому в процессе нагрева изменяется ком- 
плексный импеданс системы «индуктор - нагрева- 
емый объект», добротность и частота резонансного 
контура. Кроме того, эти величины зависят от кон- 
структивных и электрических параметров индук- 
торной системы, что особенно ярко выражено при 
нагреве ферромагнитных материалов. 

Как правило, для питания колебательного кон- 
тура используется мостовая схема резонансного 
инвертора, особенностью работы которой является 
введение системы фазовой автоподстройки часто- 
ты, параметры которой во многом определяются 
такими параметрами колебательного контура, как 
резонансная частота, добротность, а также диапа- 
зоны изменения этих величин в процессе нагрева. 
В связи с этим актуальной является задача исследо- 
вания и определения зависимостей параметров ко- 
лебательного контура от температуры, конструк- 
тивных и электрических характеристик индуктор- 
ной системы. Несмотря на существование ряда ра- 
бот [1-3], посвященных исследованию импедан- 
сных характеристик индукторной системы, вопро- 
сы расчета требуемых диапазонов подстройки ча- 
стоты и изменения добротности резонансного кон- 
тура остаются не освещенными. Решение этой 
сложной задачи является целью данной работы, в 
которой рассмотрен случай нагрева ферромагнит- 
ных заготовок при последовательном соединении 
компенсирующей емкости с индуктором. 


Математическая модель системы 
«индуктор - нагреваемый объект» 


Систему «индуктор - нагреваемый объект» 
можно представить эквивалентной схемой (рис. 1), 
которая может быть преобразована к последова- 
тельному соединению эквивалентной индуктивно- 
сти Ь э и эквивалентного активного сопротивления 
К э , импеданс которой можно найти по широко из- 
вестным формулам [1]: 

2 = Я э + ]соЬ э ; 


К 3 - к п + к д 


со 2 Ь і 


я*+ю 2 (х„+^+4) 2! 




^+со 2 (А д +А,)(А д +А, +4) 
/?;,+® 2 (С, +А, +А ) 2 


( 1 ) 


где со - рабочая частота, К и - активное сопротивле- 
ние витков индуктора, Ь ІС - индуктивность обрат- 
ного замыкания, Ц - индуктивность рассеивания, 
К л , Ь д - активное сопротивление и индуктивность 
нагреваемого объекта. 



Рис. 1. Эквивалентная схема системы «индуктор - нагревае- 
мый объект» при последовательном соединении 
компенсирующего конденсатора с индуктором 


Индуктивность обратного замыкания обусло- 
влена магнитным потоком вне индуктора, ее расчет 
затруднителен, прежде всего, из-за сложной фор- 
мы поля на этом участке магнитной цепи. Задача 
расчета Ь к для цилиндрических деталей решается 
путем введения коэффициента Нагаока к согласно 
формуле [1]: 

А, = Ѵ /с , , ( 2 ) 

1 -к 

где А 0 - индуктивность пустого индуктора: 
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4 = 


кР 2 р 0 Ж 2 

4/ 


( 3 ) 


где О и / - диаметр и длина индуктора, IV- число 
его витков, ,а 0 =1,25-1О 6 - магнитная постоянная. 

В инженерных расчетах коэффициент Нагаока 
определяется соотношением [2]: 

к(Р/1) = 0, 224 + 0, 756 • ехр( -0, 452 Р //). (4) 


Индуктивность рассеивания Ь, учитывает маг- 
нитный поток в зазоре между нагреваемой деталью 
и индуктором и рассчитывается по формуле [3]: 


_лр 0 1У 2 (Р + к)к 


( 5 ) 


где к - величина зазора. 

Активное сопротивление заготовки вычисляет- 
ся по формуле [3]: 


= И' Б \ сор(Ѳ)р 0 р(Н^Ѳ) 


( 6 ) 


К д =соЬ д . (10) 


Таким образом, импедансные характеристики 
системы «индуктор - нагреваемый объект» зависят 
от тока, протекающего через индуктор, что значи- 
тельно затрудняет расчет. 


После расчета импедансных характеристик рас- 
считывается емкость колебательного контура из 
условий резонансного режима работы: 


С = 


®;Тэо+ Л ° 


4Р 


(П) 


где Ь х , і?эо - значения Ь э и Я э при температуре 
окружающей среды. 

Попытка аналитически выразить резонансную 
частоту колебательного контура из соотношений 
(1-11) приводит к необходимости решать иррацио- 
нальное уравнение, поэтому целесообразно приме- 
нение численных методов. 


где р(Ѳ) - удельное сопротивление, 5 - площадь 
нагреваемой поверхности, р(Н 0 ,Ѳ) — магнитная 
проницаемость среды, зависящая от напряженно- 
сти поля # 0 и температуры Ѳ. 


Зависимость магнитной проницаемости боль- 
шинства углеродистых сталей от температуры 
определяется соотношением [3]: 


Р(Н 0 ,Ѳ)= 1 + 


«- Я /-1 

л V 


1 + 


уѲк -Ѳ 


Ѳ<Ѳ К 


( 7 ) 


Р = 1 , ѳ > Ѳ к 


где а=3-10 5 , /3= — 0,85, ^=1,9, 5=0,16 - коэффици- 
енты, полученные в результате компьютерной об- 
работки усредненных экспериментальных данных 
по методу наименьших квадратов [3], Ѳ к — темпера- 
тура, соответствующая точке Кюри. 

Для расчета температурной зависимости удель- 
ного сопротивления нагреваемого материала 
объекта в первом приближении можно воспользо- 
ваться классическим соотношением [3]: 

р(Ѳ) =р 0 (1 + а т Ѳ), (8) 


где р 0 - удельное сопротивление при температуре 
О °С, а т - температурный коэффициент удельного 
сопротивления. 


Напряженность магнитного поля определяется 
по формуле [3]: 


Ш 


( 9 ) 


где / 0 - амплитуда тока через индуктор. 

Индуктивное и активное сопротивления на гра- 
нице проводящей среды равны [3], поэтому индук- 
тивность и активное сопротивление нагреваемого 
объекта связаны соотношением: 


Расчет температурных зависимостей резонансной 
частоты и добротности последовательного 
колебательного контура для системы 
«индуктор - нагреваемое тело» 

При изменении температуры детали изменяют- 
ся ее электрофизические характеристики, соответ- 
ственно, меняется глубина проникновения вихре- 
вых токов в нагреваемый объект. Глубина прони- 
кновения, в свою очередь, определяет собственную 
частоту и добротность резонансного контура. Ука- 
занное обстоятельство значительно затрудняет рас- 
чет температурной зависимости резонансной ча- 
стоты и добротности, поэтому для их расчета ис- 
пользованы итерационные вычисления. На рис. 2 
показан предлагаемый алгоритм вычисления емко- 
сти компенсирующего конденсатора и зависимо- 
сти резонансной частоты и добротности от темпе- 
ратуры при постоянном токе индуктора. В блок- 
схеме алгоритма приняты следующие обозначения: 
Ѳ„ - наименьшая температура, приравниваемая к 
температуре окружающей среды, Ѳ в - наибольшая 
температура, АѲ - приращение температуры, / аж - 
номинальная частота (резонансная частота конту- 
ра при Ѳ=Ѳ Н ), / ста6 - стабилизированная амплитуда 
тока, протекающего через индуктор, / - текущее 
значение частоты, /„(Ѳ) - резонансная частота при 
заданной температуре, е - наибольшая абсолютная 
погрешность вычисления частоты, 0 - доброт- 
ность резонансного контура. 

На рис. 3 показаны зависимости резонансной 
частоты и добротности от температуры, получен- 
ные в результате вычислений по алгоритму (рис. 2). 
Зависимости построены для индуктора с числом 
витков Ж=5 и геометрическими размерами 
0=0, 1 м, /=0, 1 м, при величине зазора между витка- 
ми индуктора и нагреваемой деталью И = 5 мм, при 
/„ом=10 кГц и / ста6 = 1 кА. 
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Рис. 2. Алгоритм вычисления температурных зависимостей резонансной частоты и добротности колебательного контура при 
фиксированных значениях / па6 , Т нш , Н и И/ 



а) б) 

Рис. 3. Температурная зависимость: а) резонансной частоты ; б) добротности колебательного контура при И/=5, 0=0, 1 м, 
1=0, 1 м, Ь=5 мм, Т„ Ш =Ю кГц и І ст6 =1 кА 


Анализ полученных в результате расчета зави- таль» качественный характер зависимостей сохра- 
симостей (рис. 3) показывает, что при температуре няется. На практике температурные зависимости 
вблизи точки Кюри наблюдается скачкообразное характеризуются более плавным переходом через 
изменение частоты и добротности системы. В слу- точку Кюри, поскольку в расчетах не учитываются 
чае изменения геометрических и электрических динамические изменения теплофизических 
параметров системы «индуктор - нагреваемая де- свойств материалов. 
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Зависимости диапазонов подстройки частоты 
и изменения добротности от конструктивных 
и электрических параметров системы 
«индуктор - нагреваемое тело» 

Подстройка резонансной частоты колебатель- 
ного контура в процессе нагрева происходит авто- 
матически, поэтому для практических расчетов 
преобразователей частоты (ПЧ) важно определе- 
ние не столько зависимостей подстраиваемых па- 
раметров от температуры, сколько предельный ди- 
апазон их изменения. Эти параметры, как показы- 
вает практика, существенно зависят от конструк- 
тивных и электрических параметров системы «ин- 
дуктор - нагреваемый объект», в частности от ве- 
личины зазора между индуктором и деталью, опре- 
деляющего поток рассеяния, от тока индуктора, 
номинальной частоты в начальных условиях и т. д. 

На рис. 4 приведены графические зависимости эк- 
стремальных значений резонансной частоты и доброт- 
ности колебательного контура от величины зазора 
между индуктором и деталью, полученные в результате 
расчета по алгоритму (рис. 2) при разных значениях / т6 . 

Из приведенных результатов видно, что с ростом 
зазора диапазон подстройки частоты уменьшается, а 
границы частотного диапазона приближаются к но- 
минальной резонансной частоте. Это объясняется су- 
щественным увеличением индуктивности рассеяния, 
при котором изменение индуктивности детали Х д уже 
не оказывает заметного влияния на эквивалентную 
индуктивность всей системы. Следует отметить, что 
при проектировании индуктора стремятся макси- 
мально уменьшить зазор, т. к. это позволяет умень- 
шить требуемую мощность резонансного конденса- 
тора. Зависимость диапазона подстройки частоты от 
тока индуктора обусловлена уменьшением магнит- 
ной проницаемости стали с ростом напряженности 
поля индуктора согласно (7). Верхняя граница резо- 
нансной частоты зависит от тока индуктора, несмо- 


тря на то, что в области температур выше точки Кю- 
ри ^=1. Это объяснимо изменением резонансной ем- 
кости по условию сохранения прежней номинальной 
резонансной частоты (10 кГц). Аналогичным образом 
объясняется рост добротности с увеличением зазора: 
при Н = 0 добротность не меняется и равна единице, 
т. к. индуктивность индукторной системы почти пол- 
ностью определяется индуктивностью детали Ь д . 

Добротность может быть выражена из формул 
(1), (10) путем элементарных преобразований: 

і? д + (ь д +ьлі д +Ц+ 4 ) 

^ II 

Ь Д Ь іс 

Обычно Ь іс намного больше Ь д [2, 3], тогда 

ід гг "Г 1д ]_, 

0 = ^-^ = 1 + у ^. ( 12 ) 

ь д ь д 

Из выражений (5) и (12) следует, что при И = 0 
добротность 0= 1 и не зависит от температуры заго- 
товки в индукторе, рис. 4, б. 

Другой важной задачей является исследование 
изменения характеристик резонансного контура 
при варьировании количества витков индуктора. 
На рис. 5 приведены графические зависимости гра- 
ниц подстройки резонансной частоты колебатель- 
ного контура и добротности от количества витков 
индуктора при й=0,5 см. 

Характеристики (рис. 5, б) при малом количе- 
стве витков индуктора могут быть неточны, по- 
скольку увеличивается магнитный поток рассеива- 
ния. Можно отметить общую тенденцию сужения 
диапазона подстройки частоты с увеличением 
плотности намотки витков индуктора, что приво- 
дит к росту напряженности магнитного поля в ин- 
дукторе (9). В этом случае нагрев происходит при 
гораздо меньших значениях магнитной проница- 
емости, соответственно уменьшается и диапазон ее 
изменения. Хотя из этого не следует, что необходи- 




Рис. 4. Зависимость диапазона изменения: а) резонансной частоты; б) добротности колебательного контура от величины за- 
зора между индуктором и деталью при различных значениях / па6 и при ѴѴ=5, 0=0 , 1 м, 1=0, 1 м, Г иж =10 кГц 
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Рис. 5. Зависимость границ изменения: а) резонансной частоты ; б) добротности колебательного контура от количества витков 
индуктора при различных значениях / сга6 и при 6=0,01 м, 0=0, 7 м, 1=0, 1 м, 1 НШ =Ю кГ ц 




Рис. 6. Зависимости границ изменения: а) резонансной частоты; б) добротности колебательного контура от величины зазора 
между индуктором и деталью при разных значениях 1 НШ и при ІЛ /=5, 0=0, 1 м, 1=0, 1 м, І а3 б=1 кА 


мо стремиться к максимальному увеличению плот- 
ности намотки витков, т. к. это приводит к увеличе- 
нию напряжения на резонансном конденсаторе. 

Далее рассмотрены зависимости диапазона под- 
стройки частоты и изменения добротности от 1т и для 
разных значений начальной частоты / нж при стаби- 
лизированной амплитуде токе индуктора / инд =1 кА. 

Из приведенных графиков (рис. 6) следует отме- 
тить значительное увеличение добротности с ростом 
частоты / тм . Это объясняется тем, что при увеличе- 
нии частоты снижается глубина проникновения 
вихревых токов в заготовке, соответственно, умень- 
шается Ь д , что при прежнем значении Ц приводит, 
согласно выражению (12), к росту добротности. 

Заключение 

Предложенный алгоритм определения границ 
изменения резонансной частоты и добротности ко- 


лебательного контура позволяет определить пара- 
метры контура при изменении температуры нагре- 
ваемого объекта при фиксированных значениях ве- 
личины зазора между индуктором и объектом, чи- 
сла витков и тока индуктора, номинальной частоты. 

На основе алгоритма получены зависимости гра- 
ниц диапазонов изменения собственной частоты и 
добротности контура от конструктивных и электри- 
ческих параметров системы «индуктор - нагревае- 
мое тело». Установлено, что с ростом величины за- 
зора между индуктором и нагреваемой деталью и 
увеличением числа витков при условии постоянства 
номинальной резонансной частоты диапазон изме- 
нения частоты сужается до нуля. Максимум измене- 
ния резонансной частоты получен при нулевой ве- 
личине зазора между индуктором и нагреваемым 
объектом. Установлено значительное увеличение 
добротности колебательного контура с ростом вели- 
чины зазора и номинальной резонансной частоты. 
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